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INTRODUCTION 


Source  characterization  studies  of  earthquakes  have  generally  supported  the  view  that 
earthquakes  are  not  simple  single  point  sources  but  rather  are  composites  of  many  small  rup¬ 
tures  caused  by  source  heterogeneities  called  barriers  and  asperities  (Papageorgiou  and  Aki, 
1985a, b).  Recent  studies  of  large  earthquakes  have  modelled  them  in  terms  of  multiple  rup¬ 
tures,  each  one  delayed  by  some  amount  of  time  after  the  initial  rupture(eg.  Barker  and  Langs¬ 
ton,  1981,  and  other  references  therein).  Small  microearthquakes  have  less  dramatic  source 
multiplicity  and  directivity  effects  although  small  short-delay  multiple  ruptures  could  greatly 
affect  the  high-frequency  characteristics  of  small-earthquakes  seismograms. 

Explosion  sources,  on  the  other  hand,  are  generally  thought  of  as  being  ideal  point  explo¬ 
sive  sources,  although  many  large  nuclear  explosions  have  an  accompanying  non-isotropic  tec¬ 
tonic  component.  However,  economic  explosions,  such  as  quarry  and  mine  blasts,  are  known 
to  be  composed  of  several  delayed  blasts,  sometimes  referred  to  as  "ripple  firing".  Depending 
on  the  time  delays  between  explosions  and  the  bandwidth  of  the  recording  seismic  instruments, 
ripple-firing  effects  should  be  recognizable  in  seismic  spectra  and  may  serve  as  simulations  of 
multiple  ruptures  in  earthquakes. 

If  two  or  more  explosions  are  fired  at  nearly  the  same  location  but  with  small  time 
delays,  the  seismically  recorded  waveforms  for  each  of  the  explosions  should  be  very  similar 
although  their  amplitudes  may  vary  if  the  sizes  of  the  delayed  explosions  are  different.  The 
power  spectrum  computed  on  a  time  window  which  incorporates  all  die  delayed-explosion  sig¬ 
nals  will  exhibit  a  characteristic  modulation  or  scalloping  pattern.  The  amplitudes  of  the  scal¬ 
lops  will  depend  on  the  degree  of  correlation  of  the  different  signals  and  their  relative  ampli¬ 
tudes,  and  the  periodicity  of  the  scalloping  will  be  determined  by  the  number  and  length  of  the 
different  time  delays.  Bell  and  Alexander  (1977)  computed  spectrograms  or  sonograms  of 
simulated  multple  explosions  and  demonstrated  that  the  spectral  modulations  would  be 
observed  throughout  the  entire  wavetrain.  This  persistence  of  the  specU'al  modulations  in  all 
phases  of  the  seismograms  would  be  a  way  of  distinguishing  multiple  explosions  from  mul¬ 
tipath  interference,  such  as  modulations  in  P  spectra  do  to  P-pP  interference.  Flinn  et  al  (1973) 
demonsuated  the  application  of  homomorphic  deconvolution,  or  cepstral  analysis,  to  the  detec¬ 
tion  of  a  simulated  multiple  events.  However,  neither  of  these  studies  applied  the.se  techniques 
to  actual  delayed  seismic  events. 

The  observability  of  ripple-fire  effects  in  mine  or  quarry  bla.sts  could  have  important 
implications  for  the  problem  of  seismic  discrimination  between  economic  explosions,  small 
nuclear  explosions  and  earthquakes  recorded  at  regional  distances  with  high-frequency  seismic 
recorders.  Past  regional  discrimination  studies  have  focused  on  comparisons  of  so-called 
"discriminant.s"  usually  extracted  from  tlic  spectra  of  explosions  and  earthquake.s( Murphy  and 
Bennett,  1982;  Gupta  et  al,  1984,  Aviles  and  Lee,  1986).  The  ultimate  objective  of  these  .stu¬ 
dies  has  been  to  (ind  a  spectral  feature  which  reflects  some  intrinsic  difference  between  explo¬ 
sions  and  earthquakes  which  could  be  used  to  identify  the  source  type.  However,  it  has  proved 
difficult  to  unequivically  distinguish  between  .spectral  differences  due  to  .source  and  tho.se 
caused  by  ncar-sourcc  geology  and  propagation  path  differences.  Aviles  and  Lee  (1986)  have 
suggested  that  ripple  fire  effects  could  have  been  the  explanation  for  their  observation  of  more 
low-frequency  energy  in  mine  blasts  compaied  to  colocated  earthquakes,  although  differences 
in  near-source  geology  could  have  also  produced  these  differences.  To  date,  there  have  been  no 
studies  which  explicitly  show  the  effect  of  ripple  fire  on  high-frequency  spectra  which  could  be 
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used  to  identify  economic  explosions. 

In  this  study,  we  report  on  the  direct  observation  of  source  multiplicity  in  the  spectra  of 
explosion  data  recorded  at  regional  and  near-regional  distances.  Mine  blasts  and  earthquakes  in 
western  Norway,  recorded  at  regional  distance  from  the  new  NORESS  regional  seismic  array, 
are  compared  which  show  that  source  multiplicity  can  be  observed  in  mine-blast  spectra  but 
not  in  those  of  small  earthquakes.  Moreover  ,  we  also  show  evidence  of  source  multiplicity  in 
mine  explosions  seismograms  ,  recorded  by  NORESS,  in  Sweden  and  in  Russia  near  Leningrad 
as  well  as  in  presumed  peaceful  nuclear  explosions  (PNEs)  in  western  Russia  recorded  at  the 
NORSAR  teleseismic  array.  This  study  will  show  how  incoherent  array-stack  spectra  extending 
to  high  frequency  can  be  used  for  source  characterization  and  identification  of  source  multipli¬ 
city  in  explosions. 
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Modulated  Spectra  of  Multiple  Events 


Our  strategy  in  the  detection  of  multiple  events  is  to  identify  a  spectral  modulation  pat¬ 
tern  which  persists  in  the  spectra  of  all  phases  throughout  the  entire  seismogram.  It  would  also 
be  desirable  to  discern  the  number  and  delay  times  of  the  different  events.  However,  as  we 
show  below,  this  can  be  very  difficult  because  scalloping  patterns  for  multiple  explosions  can 
be  very  complicated,  even  for  a  small  number  of  delays.  In  this  section,  we  consider  the 
theoretically  expected  spectra  of  seismic  signals  from  multiple  explosions,  both  from  single  and 
multiple  delays,  and  discuss  the  delay  times  we  expect  from  standard  blasting  practice. 


Single  Delay.  The  equation  for  two  signals  Jt(r),  delayed  by  time  x  and  scaled  by  a,  a  real 
value,  is  given  by 

y(f)  =x(t)-»- ow(f-T).  (1) 

The  Fourier  transform  of  this  signal  is  then 

T((o)  = 


giving  the  power  spectrum 


I  y(a))  |^=  |X((0)  |^(l-haV2aco5a)x)  =  |X((o)  |^(l-»-a^)(l-t-2-^cosa>x).  (2) 

This  spectrum  exhibits  modulations  of  period  due  to  the  coscox  term.  The  logarithm  is  now 


taken  to  "whiten"  the  spectrum  (Kemerait  and  Sutton,  1982)  giving 

£.og  lK((0)|2  =  Log[lX((o)i2(l-Hx2)l  +  Lc7g(l-»-2— 5L_coswx).  (3) 

1-i-a  • 

Using  the  approximation 

Log(l+x)=x-^-l-^ - (4) 

equation  3  can  be  expanded  to  give  the  following  ; 

Log  I ^(a)) \^  =  Log[ IXfO)) +  Pcoswx  (5) 

.  cos2(ox-t-l  .  p^  .  cos3cox-t-3cos(i)x  . 

"  2  ^  2  ^  3  '  4  ^ 


where  p=2 - The  Fourier  transform  of  the  log  of  the  Fourier  transform  is  defined  as  die 

1-f-a^ 

cepstrum.  This  function  exhibits  spikes  at  x  and  multiples  of  X,  and  the  magnitudes  of  the 
spikes  diminish  with  quefrency. 

Multiple  Delays.  A  signal  with  N  delays,  with  the  Nth  delay  being  scaled  by  o.fj  (a  real 
value)  and  delayed  be  Xp,  can  be  described  by 

y(t)  =  xit)+aixit-Xi)+ajx(t-X2)+...+o.nxit-Xf^).  (6) 

The  Fourier  transform  of  this  is 
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y((0)  =  X((oXl+aic  ^“''+a2C  ■'“”*-1-  •  •  •  -f-a^^e 
and  the  log  power  spectrum  is  expressed  by 

Log  ll'((0)l2  =  Log[lX(C0)l2(Rp)l 


(7) 

(8) 


+  Log 


[  S  (a„cos(OT„+  S  a„a^cos(o(T„-x„) )  ] 

n=l  m=rt+l 


N 

where  P  =  constant  in  z.  The  approximation  in  equation  4  can  be  applied  to  expand 

n=:l 


equation  8  with  the  first  term  being 


X 


2 

1+P 


‘  N  N 

£  (a„cos©x„+  5^  a„a^cos(0(T„-T„)  ) 


(9) 


The  following  terms  contain  increasing  numbers  of  cosinisoudal  components.  The  Fourier 
transform  of  equation  8  will  produce  peaks  at  the  x„’s,  the  (x„-x„)’s  and  their  harmonics,  plus 
peaks  at  the  quefrencies  of  the  higher  order  terms  in  the  expansion.  The  amplitudes  of  these 
spikes  are  dependent  on  the  amplitudes  of  the  delayed  signals  (  the  a„ ’s). 


Expected  Delay  Times.  Standard  blasting  practice  almost  always  involves  designing  a  spatial 
and  temporal  pattern  which  can  vary  depending  on  the  appIication(  Langefors  and  Kihlstrom, 
1963;  Dick,  et  al,  1983).  There  are  three  basic  types  of  delay  series  commonly  used  (Dick  et 
al,  1983): 

Slow  or  tunnel  delays. 

Fast  or  millisecond  delays. 

Coal  mine  delays. 

Slow  delays  are  most  often  u.sed  in  tight  blasting,  such  as  in  underground  metallic  or  nonmetal- 
lic  mines,  and  are  designed  to  cause  the  rock  burden  in  one  bla.sting  hole  to  move  prior  to  tlie 
next  detonation.  Slow-delay  intervals  are  typically  on  the  order  of  0.5  to  1.0  second.  Fast 
delays  are  most  commonly  u.sed  under  less  tight  conditions,  such  as  surface  blasting.  The  mil¬ 
lisecond  delays  can  be  quite  variable,  ranging  from  25  to  50  milliseconds  for  low-period  vibra¬ 
tions  to  in  excess  of  100  milliseconds  for  longer  periods.  Coal-mine  delays  are  a  special  .series 
of  millisecond  delays  used  in  the  hazardous  conditions  of  coal  mines  and  arc  on  the  order  of 
50  to  100  milliseconds. 

Dick  et  al  (1983)  give  three  fundamental  motivations  for  using  millisecond  delays  in  rock 
bla.sting: 

1 .  To  assure  that  a  proper  free  face  is  developed  to 

enable  the  explosive  charge  to  efficiently  fragment  and  displace 

its  burden. 

2.  To  enhance  fragmentation  between  adjacent  holes. 

3.  To  reduce  the  ground  vibrations  created  by  the  blast. 

I'he  required  coarseness  of  the  fragmentation  governs  the  length  of  the  delays,  with  longer 
delays  required  for  higher  coarseness.  For  ground-vibration  reduction,  the  required  delay  inter¬ 
val  is  directly  proportional  to  the  period  of  the  strongest  vibrations.  Also,  greater  time  delays 
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are  used  for  larger  explosives. 

Another  factor  which  must  be  considered  in  blasting  is  the  geometric  arrangement  and 
number  of  blasts  in  a  sequence.  These  factors  vary  depending  on  the  requirements  of  the  given 
blasting  situation.  Moreover,  a  blast  sequence  can  be  composed  of  more  than  one  delay  lime 
or  a  complex  combination  of  different  delay  times. 

Thus,  it  is  hard  to  predict  what  kinds  of  delays  to  expect  in  mine  or  quarry  blasting 
because  they  can  vary  depending  on  the  application.  The  shortness  of  delays  that  can  be 
detected  with  spectral  analysis  is  determined  by  the  bandwidth  of  the  recording  system.  How¬ 
ever,  equations  8  and  9  show  that  the  various  periodicities  in  the  spectral  modulations  depend 
not  only  on  the  basic  delay,  x„,  but  also  the  interaction  of  delays,  If  there  are  many 

delays  of  different  lengths,  then  the  resulting  spectrum  in  equation  8  can  be  very  complex. 
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DATA  AND  ANALYSIS  PROCEDURES 


NORESS  and  NORSAR  Data.  The  data  used  for  this  study  were  recorded  at  the  NORESS  and 
NORSAR  arrays  both  located  about  100  km  north-northeast  of  Oslo,  Norway.  The  NORESS 
array  is  a  new  densely  spaced  array,  designed  for  coherent  detection  and  analysis  of  high- 
frequency  regional  phases  (Mykkeltveit  et  al,  1983)  and  consi.st.s  of  25  vertical  short-  period 
sensors  distributed  in  a  circular  area  of  3  km  aperture.  The  NORESS  array,  which  has  been 
operating  since  late  1984,  is  located  near  the  06C  subarray  of  the  older  and  larger  NORSAR 
teleseismic  array.  NORSAR  has  42  sensors  in  7  subarrays  distributed  over  a  circular  aperture 
of  about  50  km.  and  has  been  recording  data  since  1971.  The  short-period  NORSAR  data  is 
sampled  at  20  Hz  whereas  the  NORESS  array,  designed  for  the  detection  of  high-frequency 
regional  phases,  samples  data  at  40  Hz. 

Tables  1  and  2  give  the  source  parameters  of  the  earthquakes  and  explosions,  recorded  by 
NORESV  located  in  western  Norway.  The  earthquakes  are  those  reported  in  the  bulletins  of 
the  University  of  Bergen  Regional  Seismic  Network.  'I'he  locations  are  also  plotted  on  a  map 
in  F'igure  1.  The  tabulated  local  mtignitudes  were  computed  by  the  coda-  duration  method.  The 
chemical  explosions  in  Table  2  occurred  al  a  known  mine,  A/S  Titania,  and  Blasjo,  which  is  a 
dam  constsruction  site,  and  had  announced  explosive  yields  in  pounds.  The  origin  times  and 
local  magnitudes  are  those  reported  in  the  Bergen  bulletins.  These  earthquake  and  mine-blast 
locations  are  between  300  and  400  km  north-northwest  of  the  NORESS  array. 

The  locations  of  mine  blasts  in  northern  Sweden  and  western  Russia  are  given  in  Table  3 
and  are  plotted  on  the  map  in  Figure  1.  The  S2  Swedish  event  location  was  determined  by  die 
NORESS  array,  using  the  Regional  Online  Regional  Event  Processing  Package  (RONAPP), 
develojied  by  Mykkeltveit  and  Bunguni  (1984).  All  the  other  .source  parameters  were  reported 
in  the  bulletin  published  by  the  Institute  of  Seismology  at  the  University  of  Helsinki.  The  mine 
designations  are  those  assigned  to  the  different  mines  by  Helsinki  analysts  as  a  result  of  the 
application  of  their  manual  location  method.  These  events  are  located  between  700  and  1000 
kin  to  the  east  of  NORESS. 

The  locations  of  the  presumed  peaceful  nuclear  explosions  (PNE’s),  given  in  Table  4  and 
plotted  in  Figure  2,  were  taken  from  the  PDE  bulletins.  I'hese  events  occurred  before  the 
NORF.SS  array  was  operational  but  were  large  enough  to  be  recorded  by  the  NORSAR  tele.se- 
ismic  array  in  the  near-regional  distance  range  (10  to  25  degrees). 

Incoherent  Iteam.s  Short-period  waveforms  of  regional  seismic  events  can  very  complicated. 
We  show  the  waveform  characteristics  of  these  events  in  the  form  of  incoherent  beams. 
Baumgardt  (1985)  has  shown  that  incoherent  or  envelope  beams  are  useful  for  analysing  the 
overall  sh;i[ie.s  of  complicated  seismograms  recorded  at  arrays.  In  brief,  incoherent  beamform¬ 
ing  consists  of  computing  the  log-rms  amplitudes  in  adjacent,  one  second  time  windows  on 
each  channel,  starting  about  one  minute  before  the  /Vi  -wave  on.sct  time  and  extending  through 
the  seismogram  into  the  /.g  coda.  .Similar  measurements  are  al.so  made  in  two  minutes  of  noise 
background  to  Pn.  The  one  .second  log-rms  c.stimates  arc  then  averaged  over  all  the  channels 
of  NORESS  and  plotted  versus  time.  'Flic  average  noise  levels  are  estimated  over  two  minutes 
ahead  of  Pn  and  are  plotted  as  horizontal  lines.  Incoherent  beamforming  smooths  m;my  of  the 
amplitude  (Uictuations  produced  hy  UK’al  .scattering  and  noise  effects.  This  smoothing  comes 
from  the  combined  effects  of  time  averaging  over  single-channel  traces  and  stacking  the  traces 


KNSt  ().  I  IK 


6 


I'l'biuarr  1787 


SEMI-ANN.  liEP.  Source  Multiplicity  in  Explosions  SAS-  TR-ST-Ol 


over  all  channels. 

Figures  3  and  4  show  incoherent  beams  for  filtered  wavcfo'^nis  of  a  western  Norway 
earthquake  and  mine  explosion,  respectively,  recorded  at  NORESS.  The  incoherent  beams  were 
formed  after  bandpass  filtering  each  trace  in  8  filter  bands,  2. 0-4.0,  2.5-4, 5,  3.0-5.0,  4.0-6. 0, 
5.0-7.0,  6. 0-8.0,  8.0-10.0,  and  8.0-16.0  llz.,  u.sing  third-order  recursive  Butterwoith  filters. 
Each  trace  and  average-noise  e.stimate  has  been  artificially  separated  by  0,5  log-rms  units  for 
display  purposes.  AKso,  the  major  regional  phases,  Pn,  P^,  Sn,  and  Lg  are  indicated  along 
with  the  time  windows  used  for  spectral  analysis,  discussed  below. 

These  Seams  show  that  significant  energy  is  present  in  all  phases  to  the  8-16  llz  fre¬ 
quency  band.  Another  notable  feature  we  have  observed  for  western  Norway  earthfiuakes,  such 
as  the  one  in  Figure  3,  is  the  gradual  emergence  of  Pn  from  the  noise  with  increasing  fre¬ 
quency,  This  is  evident  in  F’iguie  3  where  the  Pn  doesn’t  emerge  from  the  noise  until  the  6-8 
Hz  filter.  This  tippcars  to  be  a  common  feature  of  western  Norway  and  North  Sea  earthquakes, 
as  reported  by  Mykkeltveit  (1985).  1'he  RONAPP  proces.sor  at  NORESS  consistently  underesti¬ 
mates  distances  to  these  events,  relative  to  the  Bergen  locations,  due  to  Pn  being  undetected 
and  Pg  used  instead  in  making  single-array  locations.  By  comparison,  however,  the  Pn  wave 
.stands  out  clearly  in  all  filter  bands  for  the  incoherent  beams  of  the  Bla.sjo  bla.st  in  Figure  4. 

Figure  5  shows  an  example  of  a  compressed  incoherent  beam  plot  for  one  of  the  PNE 
explosions  recorded  by  NORSAR.  In  this  ca.se,  the  incoherent  beam  was  computed  by  averag¬ 
ing  five  .second  log-rms  amplitude  e.stimatcs  over  all  available  NORSAR  suharray  channels. 
Becau.se  the  PNEs  are  located  at  teleseismic  distances  (A  >  10  deg.)  from  NORSAR,  most  of 
the  high-frequency  energy  beyond  5  Hz  has  been  attenuated.  Therefore,  only  a  0. 6-3.0  Hz  filter 
was  applied.  Also,  in  the  case  of  the  event  in  F'iguve  5,  the  first  anival  is  direct  P  rather  than 
Pn .  The  I’-coda  window,  used  in  lieu  of  the  Pn  and  Pg  windows,  is  indiciitcd.  Also,  although 
this  event  is  at  17.4  degrees  from  NORSAR,  bodi  Sn  and  Lg  are  clearly  recorded. 


Spectral  and  Cepstral  Analysis.  The  cepstrum  was  first  introduced  in  the  early  1960’s  as  a 
method  of  analyzing  com[X)sitc  signals.  It  is  defined  as  the  Fourier  uansform  of  the  logarithm 
of  the  spectrum  (  Bogart  et  al.,  1963).  'I'iiis  transform  produces  a  function  of  delay  time,  or 
"quefrency",  which  exhibits  spikes  at  the  time  differences  between  the  onset  of  the  signal  and 
any  echos  or  repeated  signals.  Cepstral  analysis  has  been  used  previously  on  seismic  data  for 
examining  multiple  events  (Flinn,  et  al.,  1973)  and  for  depth  estimation  of  events  (Cohen, 
1970;  Kernerait  and  .Sutton,  1982). 

'Fhcre  arc  a  number  of  variations  of  the  cep.strum,  including  the  complex,  power  and 
signed  cepstra.  'I'lie  complex  cep.strum  is  computed  by  taking  the  inverse  Fourier  transform  of 
the  complex  logarithm  of  the  spectrum  (Oppenheim  and  Schafer,  1975).  It  retains  all  phase 
information,  so  it  should  be  used  when  a  return  to  the  time  domain  is  necessary.  'I'he  [Xiwer 
cepstrum  is  the  magnitude  squared  of  the  I\)urier  transform  of  the  k)g  power  spectrum.  It  con¬ 
tains  no  phiise  information.  The  cosine-srjuared,  or  signed,  cepstrum  is  the  re;il  part  of  the 
Fourier  transform  rrf  the  full  symmetric  log  fK)wer  .s|>ectnmi  (Creede  :ind  Schneider,  1976).  The 
signed  cep.strum  retains  the  phase  information  from  the  frecjiicncy  domain,  so  it  may  contain 
positive  and  negative  spikes.  Ibis  |x*rmits  better  discrimination  between  spikes  caused  by 
multiple-.source  events  and  tho.se  cau.sed  by  reflection  echos.  Rcllections,  such  as  pF  or  sF,  imo- 
duce  negative  fundamental  spikes  due  to  Lloyds  mirror  effect  (Creede  and  Schneider,  1976). 
The  signed  cepstrum  was  used  in  this  study  because  no  return  to  the  time  domain  was  required 
and  detection  of  delayed  signals  which  have  the  same  CKihinty  was  desiteil. 
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In  computing  the  power  .spectra  of  the  different  phases,  an  array  stacking  procedure  was 
used  to  obtain  smooth  spectra.  The  stacking  process  cancels  out  uncorrelated  noise  and  scatter¬ 
ing  effects  in  single  channels  and  is  essentially  the  frequency-domain  equivalent  to  time- 
domain  incoherent  beamforming.  In  this  study,  array  stack  spectra  were  obtained  using  a  pro¬ 
cedure  similar  to  that  described  by  Bache  et  al  (1985)  and  Tang  and  Alexander  (1985).  The 
Fourier  power  spectrum  of  each  windowed  phase  and  the  pre-Fn  noise  is  first  computed  on 
each  array  channel.  The  window  lengths  varied  depending  on  the  duration  of  each  phase,  but  in 
general  were  between  7  and  14  seconds  for  Pn,  Pg,  and  Sn,  and  25.6  seconds  for  Lg  and 
noi.se  for  the  regional  events  recorded  at  NORESS.  In  the  case  of  the  NORSAR  recordings  of 
the  PNEs,  51.2  .second  windows  were  used  on  the  P-coda,  .9/i,and  Lg  phases  and  on  the  noise 
ahead  of  P.  These  phase  windows  are  shown  at  the  top  of  Figures  3,  4,  and  5.  The  spectra  for 
each  phase  and  noi.se  on  each  channel  were  smoothed  with  cosine  and  Hanning  windows  ,and 
the  resulting  channel  specna  were  then  averaged  across  the  arrays.  The  average  spectra  were 
then  corrected  for  the  instrument  by  spectral  division  and  for  noise  by  .subtracting  the  noise 
power  spectrum  frequency  by  frequency  from  the  power  spectrum  of  each  phase. 

Also,  for  the  NORSAR  seismograms  of  PNEs,  a  source  correction  was  made  to  the  spec¬ 
tra.  I'he  explosion  source  model  of  von  Seggern  and  Blandford  (1972)  was  used  assuming  an 
explo.sion  medium  of  granite.  The  yield  was  estimated  assuming  the  values  in  Table  4  to  be 
proportional  to  0.9  times  the  yield. 

Some  additional  processing  of  the  log  power  spectrum  was  required  before  computing  the 
cepstrum.  First,  after  correcting  the  power  spectrum  for  instrument  and  noi.se,  the  logarithm 
was  computed.  If  the  signal  spectral  density  fell  below  that  of  the  noise  at  high  frequency,  the 
last  good  value  of  the  spectrum  was  extended  out  to  the  Nyquist  frequency.  We  found  that  this 
reduced  the  ringing  in  the  cepstrum  which  occurred  becau.se  of  the  high-qucfrency  oscillations 
in  the  noise.  Second,  the  specuum  becomes  very  large  at  very  low  frequency  becau.se  of  the 
removal  of  the  instrument  response.  To  correct  this,  the  spectral  values  below  a  cutoff  fre¬ 
quency,  1  Hz  for  NORSAR  and  2  Hz  for  NORESS,  were  set  to  the  cutoff-frequency  spectral 
value.  Third,  the  mean  and  linear  trend  in  the  log  power  specuum  were  then  removed.  Finally, 
the  log  spectrum  was  reflected  about  the  Nyquist  frequency  to  produce  the  full  symmetric  log- 
power  spectrum.  This  symmetric  spectrum  ensures  that  all  of  the  information  is  in  the  real  part 
of  the  Fourier-transform  output. 
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SPECTRAL/CEPSTRAL  ANALYSIS  RESULTS 


Western  Norway  Events.  Fngure  6  shows  the  spectra  for  the  six  earthquakes  shown  on 
the  map  in  Figure  1.  Each  plot  consists  of  tlie  sjiectra  for  the  regional  phases,  f’n,  Pft ,  Sn ,  and 
Lg,  for  the  windows  shown  in  Figure  2.  The  Pg,  Sn,  and  Lg  Sj^ectra  have  been  shifted  up, 
respectively,  by  0.5,  1.0,  and  1.5  units  relative  to  the  Pn  spectra  for  di.splay  pLirpt)ses.  The 
dashed  line  on  the  bottom  of  each  is  the  pre-Pn  noise  spectrum.  Note  tliat  the  Pn  S|x;ctrum  has 
not  been  shifted  relative  to  the  noise  spectrum.  Also,  it  should  be  recalled  the  signal  S|X’ctra 
have  been  corrected  for  noise  by  spectral  subtraction.  ,Thus,  comparison  of  the  noise-corrected 
signal  spectra  and  the  noise  spectra  indicates  the  true  signal-to-noi.se  ratio  for  Pn.  Moreover, 
whenever  noise  correction  results  in  a  negative  spectral  density,  die  noise-corrected  signal  den¬ 
sity  is  set  to  die  noise  value.  Thus,  if  the  Pn  spectra  falls  on  the  noi.se  spcctum  for  several  fre¬ 
quencies,  the  signal  does  not  exceed  the  noise  level  at  those  frequencies.  However,  it  is  possi¬ 
ble  for  the  noise-corrected  signal  amplitudes  to  be  less  than  the  noise  amplitudes.  For  example, 
the  Pn  spectral  nulls  which  fall  below  the  noise  .spectrum  in  Figure  6  represent  signal-to-noi.se 
ratios  which  are  less  dian  one. 

One  common  feature  of  all  the  earthquake  Pn  spectia  in  Figure  6  is  that  there  is  no  Pn 
signal  energy  at  frequencies  less  than  4  Hz.  The  Pn  spectra  apjiear  to  be  peaked  between  8.0 
and  20.0  Hz  with  a  strong  spectral  null  somewhere  below  4.0  Hz.  This  observation  agrees  with 
the  liltered  incoherent  beam  analysis,  discussed  above  and  shown  at  the  top  of  Figure  where 
Pn  energy  only  appears  on  incoherent-beam  traces  with  filters  above  8  Hz.  The  Pg ,  Sn ,  and 
Lg  spectra  all  exhibit  energy  above  die  noi.se  levels  at  around  2.0  Hz  which  decreases  linearly 
with  frequency.  It  is  al.so  interesting  to  note  that,  with  the  exception  of  event  Q6,  which  .seems 
to  have  a  spectral  null  at  between  2  to  3  Hz,  the  Pg  and  Lg  spectra  are  almost  parallel  with 
nearly  the  same  spectral  decay  slojie,  whereas  the  Sn  spectra  appear  to  decay  at  a  slower  rate. 
In  fact,  at  frequencies  above  8  Hz,  the  Pn  and  Sn  spectra  are  very  similar  in  slojx;.  This  sug¬ 
gests  that  common  attenuation  mechanisms  are  re.S[X)nsible  for  the  Pn  ,Sn  and  Pg  ,l.g  siicctral 
decays,  and  that  the  attenuation  of  Pn  and  Sn  is  less  than  that  of  Pg  and  l.g  . 

Figure  7  shows  the  .spectra  for  die  six  western  Norway  blast.>,  3  (NMI,  NM2,  NM3)  at 
the  Blasjo  site  and  3  (NM4,NM5,  and  NM6)  at  the  Tilania  mine.  Comparing  these  spectra  with 
the  eardiquake  spectra  in  Figure  6  reveals  that  the  explosion  .s|x*clra  are  more  complicated 
The  complications  are  caused  by  modulations  or  .scalloping  in  the  sjx'ctra.  The  nulls  of  the 
modulations  occur  at  nearly  the  same  frequencies  for  the  spectra  of  different  phases,  which  is  a 
clear  indication  that  die  scalloping  results  from  source  multiplicity.  Of  course,  modulations  of 
this  kind  can  be  produced  by  multipath  pha.ses  or  depth  phases  interfering  with  the  primary 
phases.  Such  interference  may  be  expected  for  Pn,  but  not  for  Lg ,  and  the  interference  pat¬ 
terns  in  die  spectra  would  not  be  identical,  as  they  are  in  Figure  7.  Also,  the  Pn  spectra  have 
increased  signal-to-noise  ratio  at  frequencie.s  below  8  Hz  compared  with  the  eartluiuakes  in 
Figure  6.  This  result  is  also  consistent  with  the  incoherent  beam  analysis  in  Figure  4.  where 
strong  Pn  energy  is  apparent  in  all  frequency  bands. 

Figure  8  shows  the  cepslrums  for  tlie  earthquakes,  computed  by  Fourier  transloiming  the 
spectra  in  F'igure  6  using  the  methods  di.scu.ssed  in  the  |>revious  section.  Only  queliencies  out 
to  0.6  seconds  are  displayed.  'F'he  cep.strum  for  each  jilitise  has  been  shifted  in  units  of  0.05  for 
display  purposes.  'Flic  Pn  cepstium  is  unshilied.  Also,  the  parallel  dashed  lines,  labeled  for 
each  jiliase,  indicate  the  zero  levels  for  the  shifted  cepstrums  corre''poiiding  to  th.it  phase. 
Thus,  any  cepstral  value  which  falls  below  the  dashed  line  indicates  a  neg.itive  cejistral 
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amplitude.  Some  peaks  can  be  observed  in  these  cepstra,  the  strongest  ones  appearing  in  the 
low-quefrencies  of  the  Pn  pha.sc.  However,  most  of  the  cepstra  are  simple,  particularly  the  Lg 
cepstra,  which  have  nearly  zero  cepslral  values  at  all  qucfrencies. 

The  cepstrums  for  the  6  explosion  spectra  in  Figure  7  are  plotted  in  Figure  9  In  contrast 
to  the  earthquake  cepstra,  the  mine-blast  cepstra  exhibit  a  number  of  very  significant  peaks, 
indicated  by  the  superimposed  arrows.  For  a  given  event,  the  qucfrencies  of  the  major  peaks 
are  nearly  the  same,  which  is  a  manifestation  of  tlie  identical  spectral-modulation  patterns, 
pointed  out  above,  in  the  spectra  in  Figure  7. 


In  the  theoretical  discussion  above,  we  showed  that  delayed  explosions  w'ould  result  in 
modulated  sjiectra,  with  the  spacing  of  the  peaks  and  troughs  in  the  modulations  determined  by 
the  one  over  the  delay  time,  T.  For  a  single  delay,  the  power  spectrum  would  have  one  single 

modulation  of  (Equation  2).  However,  the  log  of  the  power  spectrum  would  have  higher- 

order  modulations.  The  .strongest  .secondary  modulation  would  have  periodicity  of  —X  and  its 

amplitude  would  be  negative  (Equation  5).  Thus,  the  cepstrum  for  a  single  delay  should  have  a 
positive  peak  at  a  qucfrency  of  T  followed  by  a  negative  peak  at  a  quefrency  of  2t.  There 
would  then  follow  other  .secondary  peaks  at  integer  multiples  of  T,  but  their  amplitudes  would 
decay  rapidly  with  quefrency. 


Examination  of  Figure  9  shows  that  explosions  NM2,  NM5,  and  NM6  may  have  been 
single-delay  shots.  NM5  has  main  peaks  at  about  0.04  seconds  and  a  secondary  negative  peak 
at  0.08  seconds,  or  2  time.s,  which  would  indicate  a  single  delay  of  40  milliseconds.  However, 
we  believe  that  they  may  be  processing  artifacts,  caused  by  Uie  detrending  and  rettection  of  the 
spectrum  at  20  Hz.  The  spacing  between  nulls  in  the  NM5  .spectra  in  Figure  7  is  on  the  order 
of  8  or  9  Hz  corre.s|X)nding  to  the  peak  in  the  NM5  cepstrum  in  Figure  9  at  about  O  il 
seconds,  indicated  by  the  first  arrow.  There  is  also  a  very  weak  trough,  pointed  out  by  the 
.second  arrow,  at  about  2  time.s  0.  II  second.s,  or  0.22  .seconds.  Thus,  we  would  interpret  event 
NM5  as  being  two  blasts  delayed  by  about  111  milli.seconds. 


For  NM6,  there  is  a  primary  positive  peak  at  about  0.15  and  a  smaller  negative  peak  at 
0.3  second.s,  consistent  with  a  single  delay  of  150  milluseconds.  This  seems  to  be  consistent 
with  the  NM6  spectrum  in  [-igurc  8,  which  has  at  least  two  modulations  with  nulls  separated 

by  about  fl’c  NM2  cepstra  are  the  .simplest  of  tho.se  in  Figure  9  ,  with  a  pri¬ 

mary  peak  at  about  0.08  seconds  and  the  secondary  negative  trough  at  about  0.14  seconds. 


which  suggests  a  single  delay  of  80  milliseconds.  This  corresponds  to  the  a  broad  modulation 

in  the  NM2  spectrum  in  Figuie  8  of  about  — ^ —  or  12.5  Hz. 

'  0.08 


The  cepstra  for  events  NMl,  NM.3,  and  NM4  are  more  complicated  with  more  than  one 
peak  and  trough,  indicated  by  the  solid  arrows.  From  eriuations  8  and  9  and  simulation  stu¬ 
dies,  di.scus.scd  below,  thc.se  multiple  peaks  apjicar  to  be  an  indication  that  these  explosions 
consist  of  two  or  more  combinations  of  blasts.  The  qucfrencies  of  the  primary  peaks  vary 
from  event  to  event,  but  in  general,  they  tend  to  be  between  0.08  and  0.15  seconds.  These  are 
consistent  with  multiple  explosions  with  delays  on  the  order  of  80  to  1.50  milli.seconds.  All 
these  delays  ,  including  the  single  delay  blasts,  apjicar  to  be  the  fast  or  millisccond-tyise, 
mine-blast  delays,  discussed  above. 


Swedish  and  Russian  Mine  Blasts.  I'igure  10  shows  spectra  for  the  three  mine  blasts  in 
northern  .Sweden  and  the  two  mine  blasts  near  Leningrad,  whose  locations  are  tabulated  in 
Table  3  and  plotted  in  Figure  1.  .\s  with  the  western  Norway  spectra,  the  /’g  ,  .S'li,  and  l-g 
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spectra  have  been  shifted  up  respectively  by  0.5,  1.0,  and  1.5  units  relative  to  ilie  Pn  and 
noise  spectra  for  display  purposes.  Although  these  events  are  about  twice  the  distance  from 
NORESS  as  the  western  Norway  events,  they  still  have  significant  energy  in  all  phases  out  to 
frequencies  of  16  to  18  Hz, 

The  sfiectra  of  the  Swedi.sh  events,  SM2,  SM3,  and  SMI,  are  very  complicated  with 
modulations  although  these  modulations  are  the  same  for  all  the  phases.  Since  these  modula¬ 
tions  do  not  apixrar  in  the  noise  spectrum,  they  must  be  due  to  the  signal  characteristics  and 
not  to  the  low  signal-to-noise  ratios.  Particularly  striking  in  the  SMI  and  SM3  S|)ectta  is  the 
large  broad  s|rcctral  peak  at  about  9  or  10  Hz.  SM2  also  has  broad  peak  at  about  8  Hz. 

The  spectra  of  the  Russian  blasts,  LM2  and  LMl,  are  simpler  than  the  Swedish  events 
although  spectral  modulations  can  be  easily  seen  and  are  nearly  identical  in  all  pha.scs.  How¬ 
ever,  the  spectral  periodicity  of  the  modulations  is  greater  than  that  of  the  Swedish  s[iectra  and 
more  closely  resemble  the  western  Norway  mine-bia.st  specua  in  Figure  7.  In  fact,  the  spectra 
for  tlie  BLA  mine  bla.st  NM2,  in  Figure  7,  are  very  similar  to  tliose  of  LM2.  They  both  exhibit 
two  broad  cycles  witli  nulls  .separated  by  about  12  Hz. 

The  cepstra  for  the  spectra  in  Figure  10  are  shown  in  Figure  11.  As  might  be  expected, 
the  cepstra  of  the  Swedish  mine  blasts  SM2,  SM3,  and  SMI  are  very  complicated  with  many 
peaks  and  troughs,  indicated  by  the  arrows.  SM3  and  SMI  have  strong  peaks  at  about  0.1 
seconds  quefrency,  which  probably  corresjxtnds  to  the  9  to  10  Hz  spectral  modulation  pointed 
out  above.  These  events  seem  to  consist  of  more  than  two  blasts,  although  die  basic  delay  time 
is  probably  on  the  order  of  100  milli.seconds. 

The  cepstra  for  evetit  LM2  has  a  {xrsitivc  peak  between  0.06  and  0.08  seconds,  which 
appears  to  correspond  to  the  broad  loop  with  nulls  separated  by  about  12  Hz  in  the  specba  in 
Figure  10.  There  is  also  a  trough  near  0.12  .second.s,  which  is  almost  twice  the  quefrency  of 
the  main  peak.  Therefore,  we  interpret  LM2  as  being  two  explosions  .separated  in  time  by 
about  80  milliseconds,  like  the  NM2  BLA  mine  bla,st  di.scu.s.sed  above. 

The  LMl  cepstra  contain  more  peaks  at  higher  quefrencies  than  the  LM2  cepstra.  The 
LMl  spectra  in  Figure  10  consist  of  about  two  loop.s,  the  first  having  nulls  spaced  apart  by 
about  6  Hz  and  the  second  by  about  7  or  8  Hz.  The.se  loops  should  produce  positive  cepstral 
peaks  between  0.125  and  0.17  seconds  quefrency  and  negative  peaks  at  between  0.25  and  0.34 
seconds.  As  shown  by  the  arrows  in  Figure  11,  broad  peaks  are  visable  at  these  quefrencies  in 
the  l.Ml  cepstra.  However,  sliong  peaks  are  also  vi.sable  at  quefrencies  below  0.1  seconds.  The 
origin  of  these  peaks  is  not  clear  from  the  .s[x;ctra  in  Figure  10,  since  there  do  not  appear  to  be 
broad  modulations  with  nulls  separated  by  10  Hz.  or  greater.  Processing  artifact,  jxissibly  as  a 
result  of  removing  the  linear  trend,  may  have  caused  a  humped  spectrum  of  .^0  Hz  bieadth 
which  could  have  caused  the  jxiaks  near  0.05  .seconds.  We  therefore  conclude  mat  LMl  con¬ 
sists  of  at  least  two  explosions,  .separated  by  between  125  and  170  milli.seconds. 

Western  Russia  PNEs.  'The  NORSAR  s|x'ctra  for  the  presumed  PNEs  in  western  Russia 
are  shown  in  Figure  12.  'The.se  spectra  were  computed  in  the  same  manner  as  the  NORESS 
spectra,  however,  none  of  the  spectral  plots  in  Figure  12  have  been  shifted  for  di.splay  pur- 
(X)ses.  Becau.se  these  events  are  at  distances  beyond  10  degrees,  these  s|X'ciia  do  not  contain 
the  same  high-fre(]uency  content  as  the  regional  signals  recorded  at  NORliSS.  The  NORSAR 
short-period  aniialias  filter  cuts  off  frequencies  above  5  Hz,  so  we  have  only  plotted  the  spectra 
to  5  Hz.  Clearly,  the  reduced  bandwidth  for  these  telcseismic  signals  makes  it  im[x)ssiblc  to 
resolve  millisecond  delays  in  mine  bla.sts.  Only  delays  of  200  to  300  milliseconds  or  greater 
will  be  rc.solvable  with  the  5  Hz  bandwidth  of  the  short-[X'riod  NORSAR  data. 
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Examination  of  the  spectra  in  Figure  12  reveals  modulations  like  those  observed  in  the 
mine  blasts  at  NORESS,  except  the  PNE  modulations  have  much  shorter  spectral  periodicity. 
All  the  spectra  have  a  low-frequency  spectral  null  between  between  1  and  2  liz.  The  frequency 
of  the  first  spectral  null  is  different  for  different  events,  but  the  same  in  the  spectrum  of  each 
phase  (P  coda,  Sn,  and  Lg)  for  each  event.  Spectral  modulations  are  clearest  in  the  P-coda 
spectra  which  have  the  most  high-frequency  content.  The  Sn  and  Lg  spectra  roll  off  faster 
with  frequency  and  have  less  high-frequency  content  than  the  P-coda  spectra.  Because  of  this, 
spectral  nulls  are  less  obvious  above  2  Hz  in  the  Sn  and  Lg  spectra  than  in  the  P-coda  spectra. 
Also,  as  in  the  case  of  the  NORESS  recordings  of  mining  explosions,  the  nulls  do  not  appear 
in  the  noise  spectra,  and  thus  the  signal  spectral  nulls  cannot  be  due  to  noise  contamination. 

Events  E2  and  E6  show  the  clearest  evidence  of  spectral  nulls  throughout  the  entire  fre¬ 
quency  band  in  the  spectra  of  all  three  pha.ses.  For  E2,  the  nulls  appear  at  between  1.75  to  2.0 
Hz,  about  3.5  Hz  and  near  5.0  Hz.  The  nulls  in  the  E6  spectra  occur  at  between  1.0  and  1.5 
Hz,  about  3.0  Hz,  4.0Hz,  and  near  5.0  Hz.  For  die  other  events,  after  the  first  null,  the  Sn  and 
Lg  spectra  have  one  hump  or  peak  which  is  about  the  same  as  those  in  the  P-coda  spectra. 
After  the  first  cycle,  the  modulation  pattern  in  the  Sn  and  Lg  spectra  seems  to  be  obscured  as 
the  signal-to-noise  ratio  decreases  with  increasing  frequency.  However,  in  all  cases,  there  is  an 
approximate  correspondence  in  the  location  of  the  nulls  of  the  spectral  modulations,  at  least  for 
the  first  cycle,  for  all  phases.  Thus,  the  signals  from  the.se  events  appear  to  have  been  produced 
by  more  than  one  explosion 

The  cepstra  for  these  events,  plotted  in  Figure  13.  have  peaks  and  uoughs  at  similar 
quefrencies  for  all  three  pha.ses.  However,  they  are  not  as  sharp  and  well  defined  as  those  in 
for  the  mine  blasts  recorded  at  NORESS  because  of  the  reduced  bandwidth  at  NORSAR  and 
the  fact  that  the  nulls  in  I'igure  12  are  not  as  deep  as  those  for  the  mine  blasts  recorded  at 
NORESS.  The  cepstral  peaks  below  0.4  seconds  quefrency  are  due  to  processing  artifact.  In 
particular,  the  peak  near  0.2  seconds  in  mo.st  of  the  cepstfa  is  due  to  the  reflection  of  the  spec¬ 
tra  at  5  Hz  prior  to  computing  the  cepstra. 

Event  E3  has  tlic  sharjie.st  cepstral  peaks  because  the  spectral  nulls  in  the  F3  s|iecira  in 
Figure  12  are  the  deejicst  of  all  tiic  s|iectra,  particularly  in  the  P  coda  siicctrum.  The  strongest 
P-coda  freak  for  F{3  is  at  about  1  second  quefrency.  A  strong  freak  also  appears  at  about  1 
second  in  the  Lg  cepstrum  and  a  smaller  peak  in  the  Sn  cepsinim  at  a  quefrency  slightly  less 
than  1  second.  The  slight  shift  in  the  peak  frequencies  .seems  to  be  do  to  rnnse  contamination 
in  the  Sn  and  l^g  spectra,  and  we  interpret  E3  as  being  conifxrscd  of  two  explosions  delayed 
by  about  1  second. 

The  Iv2  cepstrum  has  the  most  consi.steiit  .set  of  peaks  m  all  phases  at  ahout  0.75  .seconds 
quefrency,  as  indicated  by  the  arrow.  Event  El  has  a  set  of  somewhat  less  clear  freaks  at  a 
quefrency  of  about  0.9  seconds.  Events  1^4  and  E6  have  freaks  between  O.fr  and  0.8  seconds, 
although  are  negative.  These  troughs  may  be  due  to  noise  contamin.iiion  ami  frmecssing 
artifact. 

Invent  i;5  has  a  paiticulaily  strong  jx-.ik  in  the  P  coda  cepstrum  at  ahout  1.4  seconds, 
caused  hy  the  very  strong  spectral  modulation,  evident  m  the  l•;5  P  coda  s|x^ctium  in  Figure  12, 
with  freak  and  trough  separations  of  alunit  0  7  11/..  I  he  lirst  cycle  of  this  inodulatum  can  be 
seen  in  the  Sn  and  sfX’ctrutn,  but  then  disafrfx*ars  at  higher  Ireifuencv  due  to  noise  contami- 
n.'ition.  This  one  cycle  docs  not  apfX’ar  to  be  sufficient  to  frroduce  sliong  cepstral  peaks  in  the 
Sn  and  l.g  cefrsira  in  fugurc  13.  However,  based  on  the  strength  of  the  modulation  m  the  P- 
coda  spcctimn  and  the  fact  that  part  of  the  modulation  afrfx'ars  m  the  low  fre(|uency  ixirt  of  the 
Sn  and  l.g  sfx'ctra  m  I  .giire  12,  we  c  nclinlc  lira'  event  Efr  is  two  explosions  delayed  hy  1.4 
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SIMULATION  EXPERIMENTS 


In  order  to  obtain  a  better  understanding  of  the  spectral  characteristics  of  multiple  explo¬ 
sions,  we  simulated  multiple  explosions  using  signals  of  one  of  the  western  Norway  earth¬ 
quakes.  Event  Q1  was  used  since  its  spectrum  did  not  have  modulations  and  thus  appears  to 
be  single  event,  at  least  within  the  resolution  of  the  NORESS  bandwidth.  Thus,  if  the  signals 
from  Ql  were  produced  by  a  set  of  events  with  constant  delay  time,  the  delay  times  would 

have  to  be  less  than  — ^  Hz,  or  50  milliseconds. 

20 

Assuming  that  the  Ql  signals  are  from  a  single  explosion,  we  built  up  a  simulated  multi¬ 
ple  explosion  by  u.sing  the  Ql  signals  delayed  in  time  and  with  different  amplitude,  numeri¬ 
cally  simulating  equation  6.  'I’he  signal  was  first  i.solated  from  the  first  onset  of  the  Pn  to  well 
into  the  Lg  coda.  The  noise  ahead  of  the  Pn  onset  was  retained  and  assumed  to  be  the  noise 
ahead  of  the  multiple  explosion.  Thu.s,  the  isolated  .signal  was  added  to  the  original  seismo¬ 
gram,  but  delayed  in  time  and  multiplied  by  some  amplitude  factor,a.  This  can  be  done  more 
than  once  to  produce  a  composite  multiple-event  sequence.  We  then  computed  spectra  and 
cepstra  of  the  resultant  simulated  multiple-event  signal  as  we  did  with  the  NORESS  data. 

The  spectra  of  four  simulations  are  shown  in  Figure  14.  The  corresponding  ccpstra  for 
these  spectra  are  plotted  in  Figure  15.  Figure  14(a)  is  the  spectrum  for  two  events,  separated 
by  125  milliseconds,  and  having  the  same  amplitude.  The  resultant  .spectra  have  very  deep 

nulls  ,  separated  by  milliseconds,  or  8  Hz,  These  spectra  resemble  the  deep-null  spectrum 

of  the  Ru.ssian  mine  blast,  LMl(Figure  10),  and  the  PNE  E.l(Figure  12),  although  the  spacing 
of  the  nulls  of  E3  are  much  less  than  those  of  the  simulation.  However,  most  of  the  spectra  do 
not  exhibit  die  deep  nulls  characteristic  of  this  simulation.  The  corresponding  cepstrum  in  Fig¬ 
ure  15(a)  shows  a  clear  peak  at  0.125  seconds  quefrcncy,  as  expected.  Also,  note  the  smaller 
negative  j^xiak  at  twice  the  delay  time,  or  0.25  seconds. 

I'igure  14(b)  shows  the  effect  of  superimposing  tlrree  delayed  events,  according  to  the 
following 

.v(t)  =  a, .»•(/)  +  a2-t(f-x)  +  ayx{i-2i),  (10) 

where  (X,  =  =  1  and  x  =  0.125  seconds.  The  number  of  cycles  in  the  s|xecira  have 

been  doubled,  relative  to  the  single  delay  case,  and  the  widths  of  the  cycles  have  been  halved. 
There  are  still  two  dominant  {Xtaks,  one  at  8  Hz  and  one  at  16  Hz,  separated  by  8  Hz,  which  is 
the  reciprocal  of  the  delay  time,  0.125  .seconds.  The  ccp.stra  of  this  case,  in  Figure  15(b),  still 
have  strong  positive  peaks  at  0.125  seconds  and  an  additional  peak  at  twice  0.125  seconds,  or 
0.25  seconds. 

Figure  14(c)  shows  the  resulting  spectra  for  the  case  when  (X2  =  ~  0.5(X|  in  equation 

10,  that  is,  for  three  delayed  explosions  with  the  second  and  third  explosions  being  half  the 
amplitude  of  the  first.  The  delay  time  is  still  0.125  seconds.  Changing  the  relative  amplitudes 
of  the  different  delayed  signals  makes  the  nulls  le.ss  deep  and  the  spectral  modulations  less 
obvious.  The  cepstral  peaks  in  Figure  15  for  this  case  are  weaker  than  for  the  case  when  the 
signals  were  the  same  amplitude.  'Fhe  cejsstra  are  in  general  simpler  when  the  delayed-signal 
amplitudes  are  different,  and  the  higher  ((ucfrency  interference  |ieaks  are  much  less  apparent. 
This  simulation  re.semhles  most  of  the  cases  we  have  studied,  esi-tecially  NM5,  NM6,  and  NM2 
in  Figure  7,  I.M2  in  Figure  10,  and  all  the  I’NF-  examples,  except  F,3,  in  Figuie  12  Ituis,  it 
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appears  that  most  delayed  explosions  are  composed  of  explosions  of  different  sizes. 

F'inally,  Figure  14(d)  shows  the  case  of  30  superimposed  signals  each  delayed  by  0.125 
seconds.  Comparing  this  case  with  the  3  signal  simulation  in  Figure  14(b),  we  see  tliat  main 
peaks  at  8  and  16  Hz  begin  to  dominate  relative  to  the  secondary  lobes.  Clearly,  in  the  limit  as 
the  number  of  explosions  approaches  infinity,  the  main  peaks  become  delta  functions.  I-igure 
15(d)  shows  that  the  cepstra  for  this  ca.se  consists  of  of  the  primary  peak  at  0.125  seconds  and 
many  secondary  peaks.  None  of  the  spectra  or  cepstra  we  have  examined  so  far  resemble  this 
ca.se.  Thus,  it  appears  that  all  of  delayed  explosions  we  have  studied  consist  of,  at  most,  3  or 
4  explosions. 

It  should  be  noted  that  these  results  do  not  preclude  the  possiblility  of  there  being  shorter 
delays  than  50  milliseconds.  In  fact,  in  the  ca.se  of  the  mining  explosions,  each  explosion 
delayed  by  milliseconds  is  probably  a  salvo  of  several  charges  delayed  by  intervals  of  25  to  50 
milli.seconds  or  less  (Dick,et  al,  1983).  ITcIays  this  .short  could  not  be  resolved  by  any  of  the 
data  used  in  this  si  .  jy.  Bandwidths  on  the  order  of  80  Hz  or  greater  would  be  required  to 
resolve  delays  shorter  than  25  milliseconds. 
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DISCUSSION  AND  CONCLUSIONS 


This  study  has  shown  that  the  effects  of  ripple  firing  in  mine  blasts  and  I’NEs  can  be 
resolved  from  the  analysis  of  die  high-frequency  spectra  of  different  phases.  The  best  explana¬ 
tion  for  a  consistent  modulation  pattern  persisting  in  the  spectra  of  all  phases  in  a  seismogram 
is  that  the  signals  are  a  composite  of  two  or  more  delayed  signals.  Although  it  is  true  that 
mulitipathing  effects,  such  as  P-pP  interference,  can  also  cause  modulations,  they  would  pro¬ 
duce  different  modulations  for  different  phases.  Moreover,  it  is  unlikely  that  such  effects  would 
be  observed  in  Lg  and  Lg  -coda  spectra. 

We  have  used  cepstral  analysis  to  make  an  estimate  of  the  delay  times  between  explo¬ 
sions.  This  has  sometimes  proved  difficult  becau.>e  of  the  limited  bandwidth  available  in  the 
NORESS  spectra,  and  in  particular,  in  the  NORSAR  spectra.  However,  our  estimates  of  delay 
times  have  been  consistent  with  the  fa.st  or  millisecond  type  delays  commonly  used  in  under¬ 
ground  mines.  These  delay  times  are  consistent  with  the  known  blasting  practice  in  the 
Norwegian  mines,  where  delays  on  the  order  of  100  to  150  milliseconds  are  common  (S.  Myk- 
keltveit,  pei.sonal  communication). 

For  the  PNE.s,  our  results  indicate  much  longer  delays,  on  the  order  of  0.75  to  1.25 
second  delays.  Because  tliesc  explosions  are  much  larger  than  the  mine  blasts,  as  evidenced  by 
their  higher  body-wave  magnitudes,  they  arc  probably  being  used  to  move  much  larger 
volumes  of  rock.  'I'he  Soviet  Union  has  given  only  a  limited  amount  of  information  about  their 
PNE  program.  What  information  is  available  in  the  Soviet  literature  indicates  that  the  PNE  pro¬ 
gram  is  primarily  directed  toward  various  economic  development  applications  (Nordyke,  1983). 
For  example,  many  of  the  explosions  in  the  Pechora  River  region  and  around  die  Ural  Moun¬ 
tains,  which  are  the  regions  where  the  PNEs  used  in  this  study  occurred,  apjiear  to  have  been 
fired  near  major  river  .systems  in  order  to  redirect  their  flow  southward  through  canal  systems 
to  the  more  arid  parts  of  the  Soviet  Union  (Nordyke,  1974).  However,  recent  rcixirts  in  die 
.scientific  press  indicate  that  the  Soviet  Union  has  discontinued  these  projects  because  of  their 
likely  deleterious  impact  on  the  environment  (Dickson,  1986).  These  shots  may  have  been 
delayed  in  order  to  enhance  the  fracture  efficiency  of  the  explosions  and  al.so  to  reduce  ground 
motion  effects.  As  we  di.sciissed  earlier,  explosion  delay  times  should  increase  for  larger  explo¬ 
sions.  I  hus,  given  the  size  ot  these  explosions  and  the  magnitude  of  the  project  the  Soviets 
may  have  been  attempting,  the  0  8  to  I  25  second  delay  times  appear  to  be  reasonable 

From  visual  analysis  of  sjx'ctra  and  cc|isira,  we  have  not  been  able  to  unequivocably 
ascertain  the  number  and  sequence  of  delayed  expkisions.  The  temixiral  pattern  of  multiple- 
explosion  signals  is  determined  not  only  by  the  time  .sequencing  of  the  explosions  but  also  how 
they  are  arranged  spatially.  From  our  theoretical  and  simulation  analyses,  we  have  found  that 
more  than  two  interfering  signals  can  produce  very  com|ilex  log  spectra  and  signed  cepstra.  In 
order  to  determine  the  exac  t  pattern  of  explosions  in  our  data,  it  would  he  necessary  to  make 
simulations,  generate  theoretical  sjx’cua  and  cepstra,  and  find  those  which  niosi  closely  resem 
ble  the  observed  data.  In  order  to  infer  anything  about  the  spatial  pattern  of  cx|ilosions,  s|X’ctral 
observations  from  several  liifferent  azimuths  around  the  mine  would  be  leqmred  .Such  an 
inversion  scheme  would  be  the  sfX'ctral  equivalent  to  time  domain  studies  ot  complex  earth¬ 
quakes,  such  as  that  of  Baiker  and  l.angslon  (1 98 1). 

riicse  results  also  have  significance  for  di.sciimination  betw'cen  mine  blasis  and  eaith- 
quakes.  Although  we  have  clearly  observed  differences  between  mine  blasts  and  eaiilu(nakes  m 
western  Norway,  these  differences  are  mainly  a  consequence  of  mine  blasting  [Mactice  \ke  also 
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observed  a  spectral  difference  in  the  Pn  spectra,  where  earthquake  Pn  spectra  seem  to  be 
peaked  between  8  and  16  Hz  with  little  energy  below  4  Hz  whereas  the  explosions  have  Hatter 
spectra.  One  possible  explanation  of  this  difference  in  Pn  excitation  is  that  the  western  Nor¬ 
way  earthquakes  may  occur  at  greater  depth  in  the  crust  than  the  near-surface  mine  blasts,  and 
thus,  pPn  may  partially  cancel  Pn  energy  at  frequencies  below  4  or  5  Hz.  However,  we  do  not 
now  know  for  sure  if  this  difference  is  due  to  intrinsic  source  differences  or  to  path 
differences.  As  long  as  mine-blast  ripple  firing  is  used,  mine  blasts  can  be  distinguished  from 
earthquakes  by  looking  for  modulations  and  nulls  in  the  explosion  .spectra.  Mine  blasts  almost 
always  consist  of  delayed  explosions,  although  the  delay  times  will  vary  depending  on  what  is 
being  accomplished.  However,  we  have  not  always  observed  millisecond  delays  in  blasts  in  the 
Leningrad  region.  There  are  undoubtably  shorter  delay  explosions  used,  with  delays  less  than 
50  milli.seconds,  which  cannot  be  resolved  with  the  40  Hz  sampled  NORESS  data. 

Another  important  consideration  is  whether  chemical  mine  blasts  can  be  distinguished 
from  nuclear  explosions.  If  mine  blasting  always  uses  ripple  firing  (This  could  be  required  in 
the  negotiated  protocol  of  a  comprehensive  test  ban  treaty.),  high-frequency  data  are  available, 
and  nuclear  tests  do  not  have  delays,  then  it  can  be  done.  In  order  to  make  a  nuclear- weapons 
test  look  like  a  mine  blast,  it  would  be  incumbent  on  the  testban  treaty  violator  to  simulate  the 
ripple  fire  of  mine  blasting.  We  have  shown  in  this  study  that  PNEs  in  western  Russia  appear 
to  be  multiple  shots,  altliough  the  delays  are  greater  than  the  mine  blasts.  However,  it  may  be 
expensive  and  difficult  to  always  .simulate  very  short  delays.  Another  possibility  would  be  to 
set  off  a  string  of  chemical  explosives  and  time  a  nuclear  explosion  to  be  in  the  string 
sequence.  Such  a  blast  may  stand  out  strongly  as  compared  with  mine  blasts  which  may  be 
discernable  in  the  sj^ectra. 

Future  research  should  address  tlie  use  of  higher  frequency  data  for  the  re.solution  of 
much  shorter  delays.  Recently,  the  High  Frequency  Seismic  Element  (HFSE)  has  been  installed 
at  NORESS  (Ringdal,et  al,  1986)  which  provides  the  high  frequency  data  needed  for  such  a 
study.  The  HFSE  element  consists  of  a  Geotech  S-3  seismometer  in  a  60  m  borehole,  at  die 
NORESS  center,  and  is  sampled  at  125  Hz.  Ringdal  et  al  (1986)  have  found  that  this  element 
provides  data  which  have  good  signal-to-noise  ratio  in  Pn  and  Pg  spectra  out  to  62,5  Hz.  for 
events  at  distances  out  to  500  km  from  NORESS.  Thus,  with  the.se  data,  it  should  be  possible 
to  re.solve  explosion  delays  of  25  milli.seconds  and  less  for  the  mines  in  western  Norway. 
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FIGURE  CAPTIONS 

Figure  I.  Map  showing  locations  of  the  western  Norway  earthquakes,  the  Blasjo  (BLA) 
and  Titania  (TITA)  mines,  mines  in  northern  Sweden,  and  Leningrad-region  mines  and 
the  NORFSS  regional  array.  Locations  of  the  western  Norway  earthquakes  are  from 
tlie  bulletin  of  the  Western  Norway  Seismic  Network  of  tire  University  of  Bergen. 
BLA  and  Tll'A  locations  are  the  known  mine  locations.  'I'he  Swedisli  and  Leningrad 
mine  locations  come  from  the  bulletin  of  .seismicity  in  the  Nordic  countries  published 
by  the  Seismological  Institute  of  the  University  of  Helsinki. 

Figure  2.  Map  showing  F’DE  locations  of  the  presumed  PNEs  in  western  Russia  and  the  NOR- 
SAR  teleseistnic  array. 

Figure  3.  Filtered  incoherent  beams  of  one  of  the  BLA  mine  blasts  (NMl  in  Table  I)  in 
western  Norwtiy.  Incoherent  or  envelope  beams  are  computed  by  averaging  across  the 
NORESS  array  the  average  log-rms  amplitudes  computed  in  1  second  windows  on  each  chan¬ 
nel.  Horizontal  lines  are  the  log-rms  noise  levels  averaged  over  2  minutes  before  Pn  onset 
time.  Filters  are  third-order  Butterworlli  recursive  filters  designed  to  realize  the  bandpass  fre¬ 
quencies  indicated.  Each  plot  for  different  filters  has  been  shifted  by  1.0  log-rms  unit  for 
display  purposes.  Pn,  Pg,  Sn ,  and  Lg  windows  used  in  sp>ectral  analysis  are  shown  at  the  top. 

Figure  4.  Same  as  Figure  3  for  one  of  the  western  Norway  earthquakes  (Q6  in  Table  2). 
Diminished  Pn  amplitude  at  low  frequency  and  its  emergence  at  frequencies  above  4  Hz  are 
apparent. 

Figure  5.  Incoherent  beam  plot  for  a  we.stern  Ru.ssia  PNE  (E2  in  Table  4).  Log-rms  amplitudes 
are  computed  in  5  second  windows  on  each  channel  and  averaged  across  the  NORSAR  anay. 
'F'he  P-coda,  Sn ,  and  Lg  windows  used  in  spectral  analysis  are  indicated. 

Figure  6.  NORESS  array-average  spectra  for  each  pha.se,  Pn,  Pg,  Sn,  and  Lg,  from  the 
western  Norway  earthquakes.  The  pre-/Vi  noise  spectrum  is  plotted  with  a  dotted  line.  Each 
spectrum  has  been  corrected  for  instrument  and  pre-Pn  noise.  The  spectra  of  Pg ,  Sn ,  and  Lg 
were  respectively  shifted  by  0.5,  1.0,  and  1.5  units  relative  to  the  Pn  spectra.  The  Pn  sixjctrum 
is  not  shifted  relative  to  the  noise  spectrum. 

Figure  7.  NORESS  array-average  .s|>ectra  for  ilie  western  Norway  mine  blasts.  The  spectra 
were  computed  the  same  as  those  in  Figure  6. 

Figure  8.  Signed  cepstra  for  the  western  Norway  earthquakes,  computed  by  the  forward  Fourier 
transform  of  the  spectra  in  Figure  6.  Each  spectrum  in  Figure  6  was  detrended,  zero-meatied, 
and  rellecied  about  the  20  Hz  point  before  Fmurier  transforming.  The  Pg ,  Sn  and  Lg  cep.stra 
have  been  shifted  respectively  .05,  .Land  .15  units  relative  to  the  Pn  cepstrum.  The  dashed, 
horizontal  lines  are  the  zeio  levels  for  each  of  the  pha.se.s,  which  are  labelled  to  the  right. 

Figure  9,  Signed  cepstra  for  the  western  Norway  mine  blasts,  cotnputed  by  the  forward  Fourier 
transform  of  the  spectra  in  Figure  7  with  the  same  preprocessing  as  the  cepstra  in  Figure  8. 
Arrows  indicate  the  locations  of  peaks  which  ap|K*ar  at  the  s;ime  (picIVencics  in  all  phases  ;itid 
discussed  in  detail  in  the  text. 

Figure  10.  Spectra  for  the  noiihern  Sweden  mine  blasts  (SMI,  SM2,  and  SM3)  and  mine  blasts 
near  Leningrad  (LMl  and  LM2). 

Figure  11.  Signed  cepstra  for  the  northern  Sweden  mine  blasts  and  mine  blasts  near  Leningrad, 
computed  by  Fourier  transforming  the  spectra  in  Figure  10.  Arrows  indicate  the  quefrencies  of 
the  signilictint  cepsiial  peaks,  related  to  the  modulations  in  the  spectra  in  lagure  10,  discussed 
in  the  text. 
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Figure  12.  NORSAR  array-averaged  spectra  for  the  PNEs  listed  in  Table  4.  Spectra  are  for  P- 
coda,  Sn,  and  Lg  in  51.2  second  windows,  shown  in  Figure  5,  and  the  noise  spectium,  for  a 
51.2  second  window  before  the  P  onset,  is  plotted  with  a  dotted  line.  These  sjxctra  were  com¬ 
puted  in  the  same  manner  as  the  NORESS  spectra,  except  an  additional  source  correction  was 
also  applied,  discussed  in  tlie  text.  The  spectra  are  not  shifted  relative  to  each  other  or  the 
noise  spectra. 

Figure  13.  Signed  cepstra  for  the  PNEs  computed  by  Fourier  transforming  the  spectra  in  Figure 
12.  The  Sn  and  Lg  cepstra  have  been  shifted  up  by  0.05  and  0.1  units,  respectively,  relative  to 
the  P-coda  spectra  for  display  purposes.  The  dotted  lines  indicate  the  shifted  zero  levels  for  the 
cepstrum  of  each  phase.  Arrows  indicate  quefrencies  of  significant  peaks  related  to  the  modula¬ 
tions  in  the  spectra  in  Figure  12,  which  are  discussed  in  the  text. 

Figure  14.  Spectra  for  the  multiple-source  simulation  using  the  signals  from  the  Q1  earthquake, 
(a)  Spectra  for  composite  of  tiie  2  Ql  seismograms  added  together  but  the  second  seismogram 
delayed  by  125  milliseconds  relative  to  the  first,  (b)  Spectra  for  3  Ql  seismograms  added 
together  but  delayed  by  125  milli.seconds.  (c)  Same  as  (b),  but  with  the  second  and  third 
seismograms  reduced  in  amplitude  by  0.5  relative  to  the  first,  (d)  Spectra  for  composite  of  30 
Ql  seismograms  added  together  with  delays  of  125  milliseconds  and  with  the  same  amplitudes. 

Figure  15.  Signed  cepstra  for  the  simulated  multiple  explosions  computed  by  Fourier 
transforming  the  spectra  in  Figure  14.  The  cepstra  in  (a),  (b),  (c),  and  (d)  correspond  to  the 
s|iectra  in  Figures  14(a),  14(b),  14(c),  and  14(d),  respectively. 
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TABLE  1 

SOURCE  PARAMETERS  FOR 

WESTERN  NORWAY  EARTHQUAKES 

EVENT 

DATE 

ORIGIN  TIME 

LATITUDE 

LONGITUDE 

Ml 

(ni/d/y) 

(UTC) 

(deg  N) 

(deg  E) 

Q1 

2/05/86 

20:23:15.9 

62.66 

4.61 

2.7 

Q2 

2/05/86 

23:35:41.1 

62.74 

4.50 

2.6 

Q3 

2/06/86 

06:19:52.4 

62.90 

4.86 

2.3 

Q4 

2/13/86 

19:03:48.2 

62.61 

5.07 

2.6 

Q5 

2/13/86 

13:39:00.3 

62.40 

5.28 

2.5 

Q6 

2/16/86 

18:19:41.3 

61.69 

4.90 

2.0 

TABLE  2 

SOURCE  PARAMETERS  FOR 

SOUTHERN  NORWAY  MINE  EXPLOSIONS 

EVENT 

DATE 

ORIGIN  TIME 

SIZE 

Ml 

(m/d/y) 

(UTC) 

(TONS) 

BLASJO  MINE  EXPLOSIONS 

(  LAT.  =  59.31  deg  N  LONG.  =  6.950  deg  E  ) 

NMl 

08/05/85 

17:42:58.7 

62.9 

2.6 

NM2 

08/06/85 

17:50:07.9 

30.8 

2.4 

NM3 

10/17/85 

10:00:00.4 

32.7 

2.4 

TITAN! A  MINE  EXPLOSIONS 

(  LAT.  =  58.342  deg  N  LONG.  =  6.425  deg  E  ) 

NM4 

11/08/85 

14:18:54.6 

132.5 

2.4 

NM5 

02/14/86 

14:13:24.9 

95.7 

2.7 

NM6 

02nm6 

17:54:10.6 

16.2 

2.3 
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TABLE  3 

EPICENTERS  FOR  EVENTS  EAST  OF  NORESS 

EVENT 

MINE 

DATE 

ORIGIN  TIME 

LAT 

LONG 

Ml 

SOURCE 

(m/d/y) 

(HR-.MNtSEC) 

(deg  N) 

(deg  E) 

NORTHERN  SWEDEN  MINES 

SMI 

R1 

12/16/85 

14:44:36 

67.1 

20.6 

2.5 

Helsinki 

SM2 

01/28/86 

10:18:33 

66.0 

15.9 

2.3 

NORESS 

SM3 

R1 

02/14/86 

16:44:08 

67.1 

20.6 

2.6 

Helsinki 

LENINGRAD  MINE  EXPLOSIONS 

LMl 

VIC 

07/25/85 

17:32:46 

60.9 

29.3 

2.5 

Helsinki 

LM2 

VID 

12/25/85 

14:18:01 

60.8 

29.3 

>2.0 

Helsinki 

TABLE  4 

SOURCE  PARAMETERS  FOR 

PRESUMED  EXPLOSIONS  IN 

WESTERN  RUSSIA 

EVENT 

DATE 

ORIGIN  TIME 

LA'rrruDE 

LONGITUDE 

nth 

(ni/d/y) 

(Ul’C) 

(deg  N) 

(deg  E) 

El 

07/10/71 

16:59:59.3 

64.168 

55.183 

5.3 

E2 

10/04/71 

10:00:02.0 

61.613 

47.116 

5.1 

E3 

08/14/74 

14:59:58.3 

68.913 

75.899 

5.5 

E4 

08/29/74 

14:59:59.6 

67.233 

62.119 

5.2 

E5 

10/04/79 

15:59:57.9 

60.677 

71.501 

5.4 

F-6 

08/11/84 

18:59:57.8 

65.079 

55.287 

5.2 
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